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Using a network of seismometers and sets of optimal filters, we implemented a feed-forward control technique 
to minimize the seismic contribution to multiple intcrfcrometric degrees of freedom of the LIGO interferome- 
ters. The filters are constructed by using the Levinson-Durbin recursion relation to approximate the optimal 
Wiener filter. By reducing the RMS of the interferometer feedback signals below ~10 Hz, we have improved 
the stability and duty cycle of the joint network of gravitational wave detectors. By suppressing the large 
control forces and mirror motions, we have dramatically reduced the rate of non-Gaussian transients in the 
gravitational wave signal stream. 



I. INTRODUCTION 

The Laser Interferometer Gravitational-wave Observa- 
tory (LIGO) 1 consists of detectors located in Livingston, 
Louisiana (LLO) and Hanford, Washington (LHO). The 
goal of the experiment is to measure minute disturbances 
in the relative positions of test mass mirrors generated by 
gravitational waves, using laser interferometry. Passing 
gravitational waves are expected to produce a quadrupo- 
lar strain in the local space-time metric, resulting in a dif- 
ferential phase shift between the laser fields traveling in 
the two arms of the interferometer. Potential astrophysi- 
cal sources of gravitational waves include supernovae and 
compact binary inspirals. The expected strain from such 
sources requires a detector sensitive to displacements on 
the order of 10 -19 m/vTIz over a few km baseline, at fre- 
quencies ranging from ~*50 Hz to 7 kHz. Figure 1 shows 
a typical displacement sensitivity for the LIGO interfer- 
ometers. 

Using Fabry-Perot cavities for the interferometer arms 
as well as power recycling 2 helps improve the detector's 
sensitivity. Figure 2 shows a simplified schematic of a 
LIGO interferometer. The motion of the ground is many 
orders of magnitude above the required sensitivity, so 
the mirrors are seismically isolated, using both passive 
and active techniques. The cavities are held on optical 
resonance by suppressing external disturbances with a 
family of length and alignment feedback control loops. 
Suspending the mirrors as pendulums provides passive 
vibration isolation above the pendular resonance, which 
is arranged to be ~1 Hz. The support point of the pendu- 
lum is attached to a passive in-vacuum seismic isolation 
stack, consisting of four mass and spring layers. This pay- 
load is attached to tubes which exit the vacuum and are 
supported by an external active seismic isolation system. 
Small magnets are glued to the backs of the mirrors to 
allow for actuation via magnetic fields generated by cur- 
rents flowing in nearby wire coils. Above ~50 Hz, where 
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FIG. 1. Typical displacement sensitivity for the LIGO in- 
terferometers. LI and HI refer to the detectors at Livingston 
and Hanford, respectively. The black curve is the initial LIGO 
design sensitivity, which was not specified below 10 Hz. 

ground vibrations have been sufficiently suppressed, the 
detector is fundamentally limited by thermal noise in the 
mirrors and seismic isolation systems and by photon shot 
noise. 

There are four length degrees of freedom that need to 
be controlled. They are defined in terms of differential 
and common length changes, 
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FIG. 2. Schematic layout of a LIGO interferometer. The 
Fabry-Perot arm cavities have lengths L x ~ L y ~ 4fcm. 
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where L^, Lj,, i^, and Z y are as shown in Figure 2. An 
RF modulation scheme 3 is used to generate control sig- 
nals by relating position and angle fluctuations of the test 
mass mirrors to power fluctuations measured by photode- 
tectors at various locations around the interferometer 4 . 
The detector is brought to and held at its operating point 
by sequentially bringing cavities into resonance using dy- 
namically calculated error signals 5 , a process called lock 
acquisition. 

Precision measurement of the differential arm length 
degree of freedom, called L_ or DARM, enables the de- 
tector to find potential gravitational waves. The com- 
mon (mean) arm length, L + or CARM, is also sensitive 
to gravitational waves but is used as a reference to sta- 
bilize the frequency of the laser. When the arm lengths 
are set such that the antisymmetric port is near the dark 
fringe most of the input power is sent back towards the 
laser, such that the interferometer acts as a highly re- 
flecting mirror. By controlling the common power recy- 
cling cavity length, l + or PRC, another resonant cavity 
can be formed between the power recycling mirror and 
the Fabry-Perot Michelson, increasing the power on the 
beamsplitter and in the arms. In order to maintain the 
Michelson fringe offset required for this configuration the 
short Michelson cavity length, l_ or MICH, must also be 
controlled. 

Of the external disturbances that the control systems 
must overcome, one of the most problematic is the per- 
sistent and large 0.1 - 0.3 Hz double- frequency 'micro- 
seismic' peak, generated by wave activity in the oceans 6 . 
At the microseismic peak the amplitude is on the order 
of several pm/vHz, but the level of seismicity present at 
the observatories varies widely (by a factor of ~10) 7 . The 
relative displacement of the test masses in this band is 
of the same order as the horizontal displacements caused 
locally by the microseism. These levels of seismic noise 
present a host of problems, including: 

• Some amount of the translational ground motion 



is converted to angular motion of the mirrors due 
to cross-couplings in the mechanical dynamics of 
the seismic isolation systems. Imperfect balancing 
of the mirror actuators also converts some of the 
length control signals into angular motion. These 
misalignments create fluctuations in the power cir- 
culating throughout the cavities, driving feedback 
loops towards instability. Several of the angular 
control loops have bandwidths of only a fraction of 
a Hz. Even in the event of single arm locks the 
full interferometer may not hold resonance due to 
angular motion in the power recycling cavity. 

• Large control forces at low frequency can gener- 
ate excess noise in the signal band through non- 
linearities in the mirror actuation. Examples of 
such non-linearities are saturation of the actuator 
electronics and Barkhausen noise in the ferromag- 
nets used to actuate the mirrors. Non-Gaussian 
transients from these mechanisms pollute the grav- 
itational wave data stream as background events. 

• Some amount of light is scattered from each mir- 
ror surface, and subsequently re-enters the readout 
path after reflecting off of vibrating surfaces. When 
the relative motion of the mirror and the other sur- 
faces in the scattering path is larger than the laser 
wavelength (A ~ 1 Jim), the phase noise introduced 
by this scattered light is experienced not only at the 
frequency of the motion but as broad-band noise up 
to a cutoff frequency determined by the relative ve- 
locity 8 . 

At frequencies where the motion of the ground and the 
length fluctuations in the cavity are coherent the seismic 
noise can be actively subtracted from the cavity motion 
using a feed-forward signal. Feed-forward cancellation 
of the microseism was demonstrated at LLO during the 
commissioning of the detector, using piezoelectric actu- 
ators and filters created via system-identification 6 . Be- 
tween LIGO's third and fourth science runs (S3 and S4) 
a new external seismic isolation system, HEPI, was in- 
stalled at LLO, which used feed- forward signals to sub- 
tract local ground motion from the isolated platform mo- 
tion with a technique often called 'sensor correction' 9,1 °. 
Those techniques were used to the absolute motion of 
each individual platform. This paper will describe an 
improved feed-forward scheme which uses combinations 
of seismometers located several km's apart and Wiener 
filtering to reduce the differential motion between the 
platforms, and thereby quiet interferometric control sig- 
nals. We show results achieved during LIGO's sixth sci- 
ence run (S6) which lasted from July 2009 until October 
2010. In other work, Wiener filtering has been used for 
real-time cancellation of phase noise from vibrations in 
cavity stabilized lasers 11 , and a similar technique was 
used in a 40 meter long suspended interferometer 12 . 
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II. METHOD 

In order to optimize the noise subtraction we want to 
find the Wiener filter 13 , h, that processes the witness 
seismometer signals, w, into the best possible estimate of 
the targeted cavity length control signal, t. We find h by 
minimizing the mean square error between t and w, 
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with respect to hj, where N is the filter order (number 
of taps). This is done by setting the derivative, 
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which yields the Wiener-Hopf equations, 
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FIG. 3. Location of seismometers and vacuum chambers. 
ETMX/Y and ITMX/Y are the locations of the end and input 
test masses, respectively. RM is the power recycling mirror, 
BS the beamsplitter, MCI and MC2 the mode cleaner cham- 
bers, and STS refers to a Streckeisen STS-2 seismometer. At 
LHO the end stations were equipped with multiple single-axis 
Geotech GS-13 seismometers instead of STS-2's. 



where R is the autocorrelation matrix of the seismome- 
ter signal and p is the cross-correlation between the seis- 
mometer and the control signals. The FIR Wiener filters 
were estimated in MATLAB 14 using the Levinson-Durbin 
algorithm 15,16 , and then fit to a set of IIR coefficients to 
reduce the computational time in the real-time control 
system. The conversion from FIR to IIR was done with 
the Vectfit 17 software package. The intermediate FIR fil- 
ters were composed of 1000's of taps and processed data 
with a sample rate of 64 Hz, allowing for subtraction 
down to a few 10's of mHz. Approximately 1 hour of 
data was used to train each filter. Each test mass is 
in the vicinity of at least one seismometer (see Figure 3), 
which measures motion in three perpendicular directions, 
X, Y, and Z, where Z is vertical and the two horizontal 
directions are roughly aligned to the arms of the inter- 
ferometer. The length control is coupled most strongly 
to motion along the axes of the respective cavities, but 
cross-couplings to the other directions are non-negligible. 
Regardless of the direction of the witness signal generat- 
ing a feed-forward correction the actuation is applied in 
the direction in which the laser beam is propagating for 
that particular chamber. 

Offloading control signals to actuators located in the 
external seismic isolation systems mitigates several of the 
problems mentioned in Section I. In order to properly 
subtract the filtered witness signals, the transfer func- 
tion from our point of actuation to the cavity control 
signal must be measured and divided out. A diagram 
showing the relevant pieces of the mechanical structure 
can be seen in Figure 4. The number of mechanical com- 
ponents separating the mirrors from the ground creates 




FIG. 4. As described in Section I, the interferometer optic 
is suspended from the isolation stack. The signals from the 
ground seismometer are applied just below the stack via piezo 
(LHO) or hydraulic (LLO) actuators as indicated in the figure. 
The circles on the optic represent the four electromagnetic ac- 
tuators on the back of the mirror. The dashed line represents 
the vacuum chamber. 



a complicated transfer function with many resonant fea- 
tures. 

At LLO, the active seismic isolation system is HEPI 
(hydraulic external pre-isolator) 9,18 . A similar system 
installed at LHO, PEPI, used piezoelectric actuation. 
HEPI is an Advanced LIGO (aLIGO) system which will 
be in place at both LLO and LHO in the future. Figure 5 
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shows a schematic of the HEPI system 9 . The Wiener fil- 
ter feed-forward path was added in parallel to the exist- 
ing sensor correction path, which contained hand-tuned 
filters designed to match HEPI's position sensors to lo- 
cal seismometers on the ground. Typically this sensor 
correction reduced the differential motion sensed by the 
suspended cavities to ^ of the ground motion, in the 
microseism band. All results shown are improvements 
on top of this existing isolation. The HEPI actuators 19 
provide the ability to move its payload by ± 700 urn; the 
maximum range of PEPI was ± pm. We excited these ac- 
tuators while the detector was locked and monitored the 
mirror motion to measure the mechanical response. We 
again used Vectfit to fit these measurements. An exam- 
ple of these transfer function measurements can be seen 
in Figure G. An example Wiener filter can be seen in Fig- 
ure 7. Both the Wiener filter and mechanical response 
are only fit accurately up to few Hz. Low-pass filtering 
was applied to the feed-forward signal to prevent noise 
injection at higher frequencies. 
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III. RESULTS AND LIMITATIONS 

Before the laser beam enters interferometer, higher or- 
der spatial modes are removed by a suspended, trian- 
gular cavity called the code cleaner (MC). The length of 
the mode cleaner is controlled such that only the TEM00 
mode is resonant. The cavity is ^12 m long, with its mir- 
rors split between isolated tables in two vacuum tanks, 
MCI and MC2, as shown in Figure 3. At LLO these 
chambers are equipped with HEPI and served as the ini- 
tial testbed for the Wiener filter feed-forward implemen- 
tation. Figure 8 shows the control signal holding the 
mode cleaner on resonance with and without the Wiener 
feed-forward signal path enabled. A factor of up to ^5 
improvement was realized around the microseismic peak, 
and the RMS of the signal was reduced by a factor of ^2. 
The pre-existing sensor correction filters are left in their 
nominal state for all results shown in this paper. 

Successful implementation of Wiener filter feed- 
forward on the mode cleaner length control prompted 
commissioning of feed-forward paths for the other inter- 
fcrometric control signals, starting with the power recy- 
cling cavity control at LLO. This work was performed 
during a time of particularly high microseismic activity 
which impaired the Livingston detector's ability to re- 
main in lock. As can be seen in Figure 9(a), a significant 
reduction in control signal was achieved, with the RMS 
being reduced by more than a factor of 2. As a result 
there was also a reduction in the power fluctuations in 
the recycling cavity, which tended to drive the interfer- 
ometer control loops towards instability, and greatly im- 
proved the instrument's duty cycle. This cavity is nearly 
geometrically degenerate 20 (g- factor ^1) and therefore is 
especially sensitive to misalignments. 

Sending feed-forward signals developed to minimize 
one cavity control signal may inject noise into other 
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FIG. 5. The external isolation systems measure the motions of 
their payloads using position sensors and geophones, provid- 
ing feedback via hydraulic or piezoelectric actuators. Ground 
motion is also measured with a nearby seismometer, which is 
digitally processed and subtracted from the position sensor 
signal, a process called sensor correction. For chambers con- 
taining test mass optics this signal is then blended with the 
geophone signal into a so-called 'super-sensor' which has good 
combined sensitivity over a wider frequency band than either 
sensor individually. These 'super-sensor' signals are fed to 
the actuators. The presence of conversion matrices indicate 
a transformation between the different coordinate bases used 
by particular sensors/actuators. The use of four horizontal 
and four vertical position sensors/geophones allows for the 
sensing and control of 8 degrees of freedom: 3 translational, 3 
rotational, and 2 over-constrained. This schematic is adapted 
from Figure 3.8 of Wen 9 . 

length controls, since the cavities share mirrors. To avoid 
this problem, feed-forward paths were implemented in a 
serial fashion, such that extra motion injected into other 
degrees of freedom could be corrected by the Wiener fil- 
ters constructed for subsequent feed-forward paths. For 
example, Figure 9(b) shows the effect of filters designed 
to minimize the short Michelson control signal, calcu- 
lated on top of the existing power recycling feed-forward 
path. Above ~0.2 Hz the Michelson signal is reduced, 
with the improved isolation of the power recycling cavity 
preserved, albeit with some noise injection above ~1 Hz. 
Since the overall RMS is still reduced this noise injection 
was considered an acceptable tradeoff. 

Global feed-forward for the differential arm length con- 
trol was implemented at both observatories with similar 
results. Since the DARM degree of freedom is sensitive 
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FIG. 6. Example plant transfer function, fit using Vectfit. 
The average statistical error of points shown was ~ 2.5% in 
magnitude and ~ 1.5° in phase. 



n 

I 

*3 
s 

2f 
53 

*© 
H 
+j 
S 

o 

U 



10 



10" 



10 




—On 

- - On (rms) 

—Off 

-Off (rms) 
— Ground Motion 



1<T 10" 
Frequency [Hz] 



10 



FIG. 8. Effect of Wiener filter feed-forward on length control 
of the mode cleaner at LLO. The feed-forward signals were 
delivered to the actuators along the axis of the beam in MC2's 
HEPI system. 
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FIG. 7. Example Wiener filter, corrected for the mechanical 
plant response, also fit using Vectfit. 



to the motions of mirrors in both the end buildings as 
well as the corner station seismometers in all three loca- 
tions were incorporated into the calculation of the opti- 
mal feed- forward filter. Figures 10 and 11 show that the 
overall RMS of the control signal is reduced by a factor of 
^2.5 for each interferometer. For LLO this improvement 
in performance was balanced by a slight noise increase 
above 1 Hz. Some noise was also injected by Wiener fil- 
ter feed-forward at LHO, below 0.1 Hz, and above 3 Hz. 

The couplings of ground motion to the feedback signal 
can potentially be time dependent. Since the approach 
used here is not adaptive, this could lead to a degrada- 
tion of the filter's subtraction efficacy over time. Fig- 
ure 12 shows the performance of DARM feed-forward at 
LLO when first implemented and 8 months later. While 



pq 

1 







On 

Off 

Ground Motion 





























10" 

Frequency [Hz] 




MICH Both Paths 
MICH First Path 
PRC Both Paths 
PRC First Path 
Ground Motion 



10 

Frequency [Hz] 



FIG. 9. The upper plot shows the reduction in power recy- 
cling cavity control signal achieved using Wiener filter feed- 
forward, improving the isolation of the relative mirror motion 
to ~ °f t ne local ground motion. The lower plot shows 
the implementation of Wiener filter feed-forward on the short 
Michelson cavity and power recycling cavity simultaneously. 



the overall reduction in RMS motion was originally a fac- 
tor of ~2.5 the same feed-forward filters provided ~20% 
less isolation after 8 months, reducing the RMS motion 
by a factor of 2. There is no obstacle to retraining new 
Wiener filters to potentially recover subtraction perfor- 
mance, however if there is a change in the mechanical 
plant remeasuring the transfer function is a time con- 
suming process, due to the high accuracy required. 

As mentioned in Section I, we have observed that the 
rate of non-Gaussian events in the gravitational wave 
band increases when large forces are required to main- 
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FIG. 10. Reduction in differential arm length control from 
Wiener filter feed-forward signals at LLO. The overall RMS 
of the feedback control was reduced by a factor of almost 3, 
with some noise injection above 1 Hz. 
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FIG. 12. Performance of DARM Wiener filter feed-forward at 
the time of filter training (upper traces) and 8 months later 
(lower traces). The two sets of spectra are from times with 
significantly different seismic input. The lower plot show the 
ratio of feed-forward off/on for the two time periods. 



10 



N 

X 



es 
S 

go 
35 

"o 

t 10 
o 

U 



10 

































—On 

- - On (rms) 
—Off 
-Off (rms) 











10"' 10" 
Frequency [Hz] 



10 



FIG. 11. Reduction in differential arm length control from 
Wiener filter feed-forwad signals at LHO. The RMS is reduced 
in a simliar fashion to LLO, with some excess noise at a few 
Hz and below 100 mHz. 



tain detector lock. To measure the effect of improved 
isolation on the rate of background events, two 30 minute 
segments of data were collected for LI, with the DARM 
Wiener feed-forward first enabled and then disabled. A 
templated search, using a sine-Gaussian waveform ba- 
sis 21 , was performed and the number of triggers per sec- 
ond reported are shown in Figure 13. Only triggers with 
frequencies between 40 and 150 Hz are included, and the 
threshold SNR was set to 5. While both distributions 
show an excess of non-Gaussian events, reducing the low 
frequency control signal clearly results in a lower rate of 
signal band transients. The stationary noise floor of the 
detector did not vary significantly during the time this 



data was taken. 
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FIG. 13. Reduced low frequency control signal in the DARM 
degree of freedom also results in reduced transient rates at 
higher frequencies. This histogram shows the number of tran- 
sient events per second between 40 and 150 Hz in the LI inter- 
ferometer, over 30 minute stretches with and without global 
feed-forward. 

Fundamentally the ability to provide feed-forward sub- 
traction is limited by the coherence between the witness 
signals and the cavity control signal. A drawback par- 
ticular to the scheme detailed in this paper is the lack of 
witnesses measuring all signals which contribute to the 
cavity motion. Since the seismometers cannot distinguish 
between translation and tilt, additional sensors such as 
tilt meters could contribute additional isolation through 
more feed- forward paths 22 . There is the potential to be 
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limited by the noise floors of motion sensors, however 
at the microseismic peak the signal to noise ratio is very 
large; see Figure 14. Numerical errors can arise in the cal- 
culation of the FIR Wiener filter as well as in the fitting 
of both the IIR coefficients and the mechanical transfer 
function. These inaccuracies also have the potential to 
limit the feed-forward subtraction. 
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FIG. 14. Signal and noise levels for sensors used in the HEPI 
system. Noise floors shown for the position sensors (Kaman 
DIT-5200 IPS) and geophones (Mark Products L-4C) are the 
spectra of the horizontal over-constrained modes, i.e., the 
deformation of the beams connecting the HEPI piers. The 
seismometer (Streckeisen STS-2) noise floor was measured by 
huddle testing two sensors, i.e. positioning them adjacent to 
each other (< lm apart) such that they measure the same 
ground motion and comparing their outputs. These curves 
only represent the noise floor we are able to measure in our 
installed setup and not necessarily the noise floor of the sen- 
sors themselves. 



IV. CONCLUSIONS 

We have demonstrated the ability to use the LIGO 
detector's active seismic isolation system to reduce the 
relative motion between optics on a 4 km scale through 
seismometer feed-forward control. By sending the fil- 
tered seismometer signals to the HEPI/PEPI actuators 
we were able to quiet the relative motion of the mirrors, 
offloading control signal from the test mass actuators to 
the seismic isolation system. During times of elevated 
microseismic activity this allowed for unprecedented in- 
terferometer stability. The rate of audio band transients 
is significantly reduced by having less motion and less 
applied force. This method will be directly applicable to 



the second generation gravitational-wave detectors. 
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